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To understand the range of close-to-classical critical behavior seen in various electrolytes, gener- 
alized Debye-Hiickel theories (that yield density correlation functions) are applied to the restricted 
primitive model of equisized hard spheres. The results yield a Landau-Ginzburg free-energy func- 
tional for which the Ginzburg criterion can be explicitly evaluated. The predicted scale of crossover 
from classical to Ising character is found to be similar in magnitude to that derived for simple fluids 
in comparable fashion. The consequences in relation to experiments are discussed briefly. 
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How can one understand the fact that some elec- 
trolyte solutions display classical (or van der Waals) 
critical behavior down to deviations from criticality of 
\t\ = \T - Tc\/Tc - lO-'^ or less while others ex- 
hibit purely Ising- type criticality p] or, in some cases, 
crossover to Ising character at scales ~ lO^^'^-lO^^'^ 
The system triethylhexylammonium triethyhexyl- 
boridc (N2226B2226) in diphenyl ether so far reveals no 
hint of Ising character and is also the one that appears 
to approximate most closely the simplest sensible theo- 
retical model of an ionic system, namely, the restricted 
primitive model (RPM), consisting of = N+ + N^ hard 
spheres of diameter a with carrying charges +q and 
7V_(= A''-!-) charges —q, in a medium of dielectric con- 
stant D. Real solutions and molten salts deviate from 
the RPM in many ways: soft cores, differently sized ions, 
non-additive ionic diameters, long-range van der Waals 
interactions, short-range attractions, ionic polarizability, 
specific ionic chemical bonding, the molecular structure 
of the solvent, etc. |^|-^|. At least some of these fea- 
tures must be responsible for the observed variation of 
tx, from tx — 1 down to, perhaps, tx ~ 10~^. Neverthe- 
less, the most appealing scenario is that the RPM itself, 
as an "extremal model," displays purely classical critical 
behavior or, failing that, has a very small value of tx that 
is increased by more realistic interactions. (To match the 
results for N2226B2226 one needs tx ^ 10"^ g.) Here 
we address this scenario. Sadly, perhaps, our analysis, 
which encompasses all the leading physical effects [Q, 
does not support this attractive picture! Instead, it pro- 
vides grounds for believing that tx for the RPM is of 
comparable magnitude to that for simple molecular flu- 
ids or liquid mixtures that exhibit little if any non-Ising 
behavior. 

To introduce our approach note first that all cur- 
rently available theories predicting ionic criticality (for 
dimensionalities c? > 2 [gj) yield classical behavior be- 
cause, at heart, they are of mean-field character l^j-^j. 
Recall also that the standard approximate integral equa- 
tions such as the HNC, YBG, MSA, etc., yield no ac- 
count of the critical region or else fail to predict divergent 
critical-point density fluctuations [§j5[|^. Second, direct 
simulations of the RPM are far from being able 



to distinguish between classical and Ising criticality. In 
principle, a renormalization group (RG) treatment of the 
fluctuations could reveal the true nature of Coulombic 
criticality; but, in practice, that requires an appropriate 
LGW effective Hamiltonian which has not been available 
Furthermore, quantitative aspects become important 
if, as seems likely [5, 6(a)], a scale tx is present. 

However, as emphasized previously p|, a sufficiently 
good mean-field free energy can provide a foundation for 
an LGW Hamiltonian; then one may estimate the do- 
main of validity of classical critical behavior by using the 
Ginzburg criterion ||Tl|] which, indeed, is implied by RG 
theory If m(r) = [p(r) — pJ/Pc is the order pa- 

rameter (the overall ionic density being p = N/V) then, 
omitting the ordering field ft-, the expected LGW form in 
d spatial dimensions is 



n/kBT = a-^ j rfr[-/(m) + ^bli-^mf 
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with the spatially uniform reduced free energy density 
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The further gradient terms neglected in (|I]) and the cor- 
rections in (^ are not needed for a leading order de- 
scription of criticality: in principle, however, large or 
anomalous values of these terms could prove quantita- 
tively significant. 

To quantify the Ginzburg criterion we examine the 
fluctuations of the order parameter, normalized by the 
spontaneous order, toq {t < 0), in a d-sphere, S, of ra- 
dius set by the correlation length ^(T), that is 



dr 



(m(r)TO(O)) - (m)2 J^drGpp{r) 



(3) 



where Gpp{r) — (p(r)p(O)) — (p)^ is the density-density 
correlation function, while |S|ci=3 = ^t^^^- Since Gpp{r) 
decays fast on the scale ^, it is convenient to extend the 
integral to 00 and accept 

3x{T)/47Tp,ml{T)e{T), (t < 0) (4) 

where X = J drGpp{r)/p is the reduced susceptibility. 
Treating (|l|) and (^) simply as a mean-fleld free-energy 
functional, F[p{r),T], yields the asymptotic relations 
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for t — > 0— at p = Pc- We may now set Q = I to obtain 
an explicit estimate for the crossover scale t^, namely, 



to = (9ul/87^2c2)(a/62)^ 



for d 



(6) 



below which fluctuations dominate and the mean-field 
theory loses validity: Ising behavior should be exhibited 
for |<| < to- 

Now, judging by concordance with the simulation es- 
timates of Tc and pc [10, 8(b)], the most successful avail- 
able theory for the RPM critical region ^ is based on 
the Debye-Hiickel (DH) analysis |13|, supplemented by 
(i) Bjerrum (+, — ) ion pairing (Bj), (ii) solvation of the 
dipolar pairs in the fluid of free ions (DI), and (iii) hard- 
core repulsions (HC) P,p^. From the free energy F{p, T) 
of the DHBjDIHC theory (which is subject to minor nu- 
merical variants [4, 8(b)]) one can, by careful numerics 
[ p5[ , extract, in addition to Pc and T^, ci and U4 in (|^); 
note that /13 plays no role in (^. 

But the Ginzburg analysis demands also the coefflcient 
&2 of the gradient-squared term in ( |l|) that sets the am- 
plitude of the correlation length [scc^)]. Indeed, cor- 
responds to the range of the effective density-density at- 
tractions in the RPM: these are embodied in the density 
correlation functions and, more explicitly, in the wave- 
vector dependent susceptibility 



X(k) = 
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The last relation (for t 0+, p = Pc) follows by identi- 
fying, as before, the free energy for a nonuniform system 
with H Q. In the past, DH theory has been regarded 
as throwing light only on the charge correlation function, 
Ggq, while remaining silent on Gpp. However, in |^ we 
have shown how DH theory and its necessary DHBjDIHC 
extensions |^ can be generalized in a rather straight- 
forward way to yield a classical functional of p(r) and 
thence an explicit mean-field expression for x(l^)- Fur- 
thermore, the unexpected divergence of ^(p, T) predicted 
when p ^ turns out to be universal and exact\ By this 
route, therefore, the critical coefficient 62 is revealed and 
the crossover scale to can be explicitly calculated JlSt . 

Before describing our results for tQ, however, we note 
that Leote de Carvalho and Evans [|l6| have recently 
demonstrated the strategy set out above ||] by ap- 
pealing to the generalized mean-spherical approximation 
(GMSA). This ingenious, OZ-based |jl^ approximate in- 
tegral equation repairs the simple MSA (for which 
the density fluctuations remain bounded) by adding to 
the direct correlation functions a term with parameters 
which are adjusted to satisfy various desirable sum rules; 
thence ^ diverges at criticality and 62 can be estimated. 
Unfortunately, however, the GMSA exhibits some seri- 
ous defects: (a) the correlation length ^(p, T) varies non- 
universally and quite incorrectly when p [H , thence 



casting doubt on the plausibility of the results for p ~ 
b) the predicted value of Tc is significantly too high 
(c) no account is taken of Bjerrum pairing; and 
(d), apparently as a result of this, the GMSA free energy 
violates Gillan's upper bound throughout the criti- 
cal region (while DHBjDIHC theories satisfy it). Conse- 
quently, although the GMSA values for tc provide 
an interesting benchmark, they are surely not adequate 
for the purpose at hand. 

For reference we start with pure DH theory [^M] which 
yields C2 = l/647r, = l/30727r while feS/a^ = 
(1 -h f ln2 - 61n|)/647r ~ (0.052)^^. Via (|) these 
yield tc — 12.90. This number is large and certainly not 
suggestive of any significant regime of classical behavior; 
but the derivation of (|^) entailed various essentially arbi- 
trary numerical assignments. For calibration, therefore, 
it is essential to calculate tc by a comparable procedure 
for a simple-fluid model that one can be confident ex- 
hibits typical Ising behavior. 

To that end we start with the functional generalization 
of the Mayer expansion for a single-component fluid with 
a short-range pair potential u{r), namely, 

FS'^[p{v)]/kBT ^- j drp(r)[lnA3p(r) - 1] 
- i / drrfr'p(r)/(u(r - r'))p(r') + 0{p% (8) 



where f{x) = exp(— x/fesT") ~ 1 and A = h/y/2TrmkBT 
[8(b)]. This second- virial level suffices to describe the at- 
tractions driving criticality, for which we take a square 
well of range Xa and depth e. For the repulsions we 
adopt hard cores of diameter a: to treat these we follow 
our RPM approach and approximate the 0{p'^) terms 
by a local expression of free-volume (FV) or Carnahan- 
Starfing (CS) form |,|]. Note, however, that in treating 
both this SqWHC model and the RPM, the CS expres- 
sion is not obviously preferable since the attractive in- 
teractions (direct or effective) necessarily enter the true 
higher-order virial coefficients and act to soften the hard 
core effects. 

In the FV approximation the critical parameters are 
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with B* = l/a^pmax- For the range A = 1.4-1.7 that 
reasonably models real simple fluids |2^, Tc depends 
strongly on A. However, fc^Tc = 1.48£ should describe 
the corresponding van der Waals/classical theory quite 
well thus for the assignments B* = |7r and |\/3, 
which correspond to the exact hard-core second virial 
(2V) coefficient and bcc close packing, respectively 
we choose A ~ 1.65 and 1.43. The LGW parameters in 
(|) are found to be U4 = 3/16B*, 



C2 = {57B* + 16TTX^)s/108kBTcB* 



(10) 
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Note that in the infinite range Kac-Baker hmit, A ^ oo 
with Tc fixed, one has 62 — > oo and correctly finds to 0. 

Using the CS hard-core form one can derive a quin- 
tic equation for the critical density HJ] which yields 
Pctt^ ~ 0.249129. Normalizing to ksTc/e = 1.48 as above 
leads to A ~ 1.55. 

The third column of Table I lists the values of to for the 
SqWHC model found using the various hard-core approx- 
imations, with A chosen as indicated, and, for reference, 
with A = 1.50. For completeness, an RPA treatment ||l^] 
is included. The crucial LGW coefficients are also given. 
The last column presents the correlation length ampli- 
tude, defined via i{pc,T) w /t^^^ as i ^ 0+ 
These estimates all lie quite close to 0.41a. On the other 
hand, the values of to prove very sensitive to the approx- 
imations, ranging from 0.33 to 2.4 (even discounting the 
HC/bcc value). Since to ^ b\ [see (||)] a strong depen- 
dence on A is not surprising; but one might have hoped 
for better agreement among the approximate methods. 
Nonetheless, we may conclude that to = 10*°'^ will char- 
acterize fluids that display only Ising behavior (or 'im- 
mediate' crossover). 

We may now assess the data given in Table II for the 



RPM. In addition to ta and /a (in the second and last 
columns) it is instructive to examine the estimates for 
T* = kBTcDa/q^ and p* = PcCi^: these provide a mea- 
sure of the merit of the various approximations relative 
to the simulation data |Tc| ] which may be summarized by: 
IO^t; = 5.2-5.6, lOV* = 2.3-3.5 [8(b)]. Although the 
LGW coefficients here are factors of 3-100 smaller than 
for the SqWHC model, the ratio £,q /a remains of order 
unity and is again fairly insensitive to the approxima- 
tions: those yielding (T*,p*) in the simulation range sug- 
gest ~ 0.80. For reference, the table also lists the crit- 
ical values of Z = p/ pksT and the inverse Debye length 
K — {'i-Kq^ pi / DkuTY/'^ , where pi — p++ p^ = p — 2p2 is 
the density of free-ions while p2 is that for the ion pairs 
[H,||. The ratio {P2/ P*)c measures the degree of pairing 
in the critical region: it is quite significant [^,|l8|. 

The most striking feature of Table II, however, is the 
evidence that tc for the RPM lies in the range 10° '^ to 
10^'* and so is significantly greater than the value of Iq 
for the hard-core-square-well model! From this perspec- 
tive, the RPM should not have an unduly small region of 
Ising-like character but rather one of the same order, or 
even larger, than in simple fiuids, a conclusion certainly 
at variance with the most natural interpretation of the 
experimental evidence . 



TABLE I. Ginzburg crossover scale, tc, and critical parameters predicted for a hard-core square- well fluid (range Aa) [15]. 
See text for hard-core (HC) approximations and Ref. [16] for RPA. 



HC 


A 


ta 






C2 


U4, 


62/a 




bee 


1.433 


0.097 


1.48 


0.433o 


2.0235 


0.2436 


0.5506 


0.387 


2V 


1.65i 


1.62i 


1.48 


0.1592 


0.4926 


0.0895 


0.3234 


0.46i 




1.50 


2.418 


1.13i 


0.1592 


0.5395 


0.0895 


0.298o 


0.406 


CS 


1.553 


0.284 


1.48 


0.249i 


0.9478 


0.113i 


0.419i 


0.43o 




1.50 


0.33o 


1.334 


0.249i 


0.9792 


0.113i 


0.4066 


O.4I1 


RPA 


1.50 


1.57 


1.267 


0.2457 


0.6735 


0.1125 


0.3332 


0.406 



TABLE IL Ginzburg crossover scale, to, and critical parameters predicted for the RPM at various levels of approximations 
[4,15]: DH, pure Debye-Huckel [4,13]; GMSA [16]; DHBj, with naive ion-pairing [5,8]; +DI with dipole-ionic fluid solvation 
(and ai = a, a2 = 1.16198a [8]); hard-core treatments /bcc/2V/CS, see text and [8]; DI', with a new charging process [23]. 



Approx. 


ta 




10 


Kca 




lOZc 


10^02 


10^ M4 


62 /a 




DH 


12.9o 


6.25 


0.497 


1 





O.9O4 


0.497 


O.IO4 


O.O5I7 


0.733 


GMSA 


1.08 


7.858 


1.448 


1.522 







1.616 


0.33g 


0.0953 


0.75o 


DHBj 


12.9o 


6.25 


4.5I7 


1 


2.OI0 


4.549 


4I.OI4 


704.6 


0.4694 


0.733 


+DI 


5.36 


5.74o 


2.778 


1.123 


l.lOi 


2.236 


2.295 


1.53o 


0.1138 


0.75i 


+DI/bcc 


10.67 


5.542 


2.594 


1.029 


I.O64 


2.484 


2.223 


2.249 


0.1159 


0.777 


-I-DI/2V 


23.25 


5.227 


2.443 


0.923 


1.044 


2.823 


2.178 


3.7I7 


0.1208 


O.8I9 


-KDI/CS 


21.52 


5.249 


2.454 


0.93i 


1.046 


2.798 


2.I85 


3.548 


O.I2O4 


O.8I5 


+DI' 


1.23 


5.969 


2.38i 


I.O69 


O.9I9 


1.680 


1.77o 


7.98i 


0.1223 


O.9I9 


-hDI'/bcc 


2.25 


5.784 


2.145 


0.978 


0.852 


1.874 


1.653 


10.249 


O.I2I5 


0.945 


-hDl72V 


4.3i 


5.506 


I.9I9 


0.877 


0.79i 


2.137 


1.543 


I4.O80 


0.1226 


0.987 
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Although the DHBjDIHC theories account well for the 
leading physical effects near criticality and all leading 
terms have been included in Hlgw, it is possible that 
to, as calculated here, is a deceptive measure of the true 
RPM crossover scale, t^- Perhaps the ion-dipolc and 
dipole-dipole interactions, neglected as of order p^, play 
a special role ; this is being studied within a DH-style 
approach ||2^ . In principle, strong asymmetric terms, 
such as /13 and U5 in (||), could, under the full nonlin- 
ear RG flow, invalidate the perturbative Ginzburg anal- 
ysis. Higher-order gradient terms in Hlgw, especially 
if negative as arguments of Nabutovskii et al. |24| sug- 
gest, might, instead, bring RPM criticality within the 
crossover domain of some multicritical point ||] . The fact 
that the crossovers seen experimentally are much sharper 
than standard (taking place in a decade or less [pp6|) 
supports this view. To justify such a scenario, however, 
seems to demand a more sophisticated and quantitative 
RG analysis than normally feasible. 

Conversely, if the RPM itself does exhibit no signif- 
icant crossover from classical behavior, as our analy- 
sis indicates, the anomalous experimental results 
must be ascribed to one or more of the features lack- 
ing in the RPM that were listed initially 1^-01 . Some 
of these, like the presence of short-range van der Waals 
or solvent-mediated attractions, can and will be incorpo- 
rated in our formalism although, in truth, it is hard to see 
how they will significantly alter the to values. Indeed, 
the experimental trends seem, as mentioned, to indicate 
that ix always increases when the Coulombic forces have 
to compete with solvophobic effects. Unfortunately, one 
must also allow that particular impurities might distort 
the data in unexpected ways: one can imagine selective 
binding leading to 'big' dipoles, or long ionic 'rods.' The 
discovery, discussed in [Q , that the dimensionless correla- 
tion length parameter ^qpI^^ fits our calculations rather 
well when the Ising-fitted amplitudes are used may point 
in this direction. 

Interactions with Professor R. Evans and Dr. R. J. F. 
Leote de Carvalho have been appreciated. NSF Grant 
No. CHE 93-11729 supported our research. 
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